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Abstract 



The Lorentz Integral Transform (LIT) method is reformulated via the Lanc- 
zos algorithm both for inclusive and exclusive reactions. The new technique 
is tested for the total photoabsorption cross section of 3 H and 4 He. Due to 
the rapid convergence of the algorithm one has a decrease in cpu time by 
two orders of magnitude, but at the same time an excellent agreement with 
the results of a conventional LIT calculation. The present work opens up the 
possibility of ab initio calculations for inclusive and exclusive processes for 
A>6 with inclusion of complete final state interactions. 
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I. INTRODUCTION 



The study of reaction cross sections is an important tool to reveal the dynamics of 
particle systems. For systems with a small number of particles one aims at microscopic 
calculations trying to take into account all relevant degrees of freedom of the considered 
process. However, calculations via the classical approach, where one employs initial and 
continuum state wave functions, are very difficult to perform for reactions at energies beyond 
a possible three-body breakup threshold. On the other hand the calculation can be carried 
out in an alternative way, namely via the Lorentz Integral Transform (LIT) method [|TJ. 
This method allows ab initio calculations which take into account final state interactions 
correctly, without explicit knowledge of the complicated continuum wave functions. Various 
applications of this method for reactions of nuclei up to four nucleons can be found in the 
literature (see e.g. 

The LIT is expressed in terms of Lorentz states which have bound state like asymptotic 
boundary conditions that are much simpler to treat than the boundary conditions for a 
continuum state and thus \l/ can be calculated with similar methods as a bound state wave 
function. Common to many bound state methods is an expansion of the wave function 
on some basis set. In fact in the above mentioned LIT calculations an expansion on the 
hyperspherical harmonics basis has been employed for the reactions of nuclei with A > 3 
(with the exception of Ref. @). As the Hilbert space increases with the number of particles 
the calculations face larger and larger hamiltonian matrices, which have to be inverted in 
order to solve the equations associated with the LIT method. Such an inversion is time 
consuming and limits the number of possible basis states used in the calculation. The 
aim of the present paper is to point out that the inversion can be avoided. Inspired by 



recent condensed matter calculations |8|-[T0| we have adapted the Lanczos method for the 
calculation of the LIT. In this way we open the possibility to extend the LIT method to 
problems, where a rather large number of basis functions is needed, e.g., calculations of 
reactions of the a particle with realistic interactions or calculations of reactions of systems 
with more than four particles. 

The outline of the paper is as follows. In Sections II/III we review the LIT method for 
inclusive/exclusive reactions briefly and discuss the corresponding implementation of the 
Lanczos technique. Various checks of the calculation are discussed in Section IV, which also 
contains a short summary at the end. 



II. INCLUSIVE PROCESSES 



As pointed out for the case of the Stieltjes transform |TT| the integral transform method 
can be used for the calculation of a variety of reaction cross sections. In this and in the 
following section we restrict the discussion to inclusive and exclusive cross sections due to 
an external probe, though the approach can also be applied to other reaction cross sections, 
e.g., to the elastic scattering of particle systems. 

In order to determine inclusive cross sections due to external probes one has to calculate 
various dynamical functions 

F kV (u,q) = [ # / (Vo|^(g)|V/)(V/|Ofc'(?)IV'o)^/-.^-w), (1) 
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where u and q are energy and momentum transfer; ijj and ipf denote ground and final state 
wave functions of the considered system with energies E and Ef, respectively, while the 
various O k are excitation operators inducing different types of transitions. 

For few-body reactions with A > 2 one faces very often the problem that F k y(u, q) cannot 
be calculated exactly, since the microscopic calculation of if)f is too complicated. However, 
via the LIT approach the problem can be reformulated in such a way that the knowledge of 
ipf is not necessary [ffj. To this end the LIT of the response function is introduced 

L kkl {o)=j^ Fkk f 2 { 2 , (2) 
J (uj - a R y + of 

where 

o = o R + 207 o R , 07 > 0. (3) 

In (|^), as well as in the following equations, we suppress the dependence on q of L kk '(o,q) 
and F kk ,(u,q). 

The LIT method proceeds in two steps. First L kk i(o) is evaluated in a direct way, which 
does not require the knowledge of F kk i, and then in a second step the dynamical function 
is obtained from an inversion of the LIT. The direct calculation of L kk >(o) is outlined in 
the following. Starting from the definition of the LIT in (|2]), substituting for F kk i(uj) the 
expression in (|l]) and using the completeness relation of the Hamiltonian eigenstates, 

Jdf\^ f )(^ f \ = l, (4) 

one finds 

L w (a) = M)|6 U - 1 — > 1 — 4#o>. (5) 

ri — Eq — (Tr — to 1 H — E — Or + ioj 

It is evident that the solutions of the equations 

(jy _ e - o R + iO!)\^) = O#o) ; I = k, k' (6) 

lead directly to the searched transform: 

L kk '{a) = (7) 

Physical solutions for the Lorentz states of (0) have asymptotic boundary conditions 
like a bound state. Moreover the solutions are unique, since due to the hermiticity of H the 
homogeneous equation 

(H-E o -<T R + ia I )$ l ) = (8) 

has only the trivial solution \ipi) = 0. 

The inversion of the calculated L kk >(o) leads to the dynamical functions F kk r(u>). For 
details on the inversion procedure see, e.g., Ref. 0. The LIT method leads to reliable results 
as shown in test calculations for two- and three-body systems |],0]. 
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If one solves (§) expanding \tpi) on a set of basis functions \fj) one has to invert the 
corresponding hamiltonian matrix H^^. For a large number of basis functions this is a time 
consuming calculation, moreover, for dense matrices this task rapidly becomes prohibitive. 
As we will show, the Lanczos method can help to overcome this problem. 



Basically the Lanczos method |12[ is used to tridiagonalize matrices. If one is interested 
in calculating the ground state energy of an hamiltonian matrix H^^ one has first to choose 
a starter, a vector \<f>o), which must have a non-zero overlap with the ground state |^o)- O ne 
can then build the Lanczos orthonormal basis {\4>i),i = 0, . . . ,n} by applying the Lanczos 
algorithm recursively 

b i+1 \<f> i+1 ) = H\<pi) - Oi\(j)i) - , (9) 

where 

|0_i)=O, (hlfa) = 5ij. (10) 
The Lanczos coefficients a, and hi are defined as 

Oi = (&|ir|&), bi =\\ bt\<f>i) \\ , b = 0. (11) 

If one applies M Lanczos steps to an M x M matrix, then the original matrix is reduced 
into a tridiagonal form. However, the power of the Lanczos algorithm lays in the fact that 
even after m « M steps the eigenvalues of the m x m tridiagonal Lanczos matrix are 
good approximations of the extreme eigenvalues of H. These eigenvalues converge rapidly 
with increasing number of Lanczos steps. Working with finite precision, it is necessary to 
recall that at each step there are round-off errors so that after a certain number of steps the 
Lanczos vectors loose their orthogonality. In order to obtain accurate results this loss has 
to remain small, otherwise one has to apply an additional reorthogonalization procedure. 

In the following we show how the Lanczos algorithm can be applied to the LIT method 
for inclusive processes. To this end we rewrite the LIT in the following form 

-L w {a) = —1m\(ik>\dt — !— grOfc/^o)) • (12) 

A similar relation connects the dynamical functions F kk / to the Green's function 

F kk ,(uj) = --Im\\imG(LU + ir ] + E )} ; G(z) = (MOl—^O^o) , (13) 

provided that z = u + irj is replaced by ctr + zcrj. This is not surprising since for 07 — > the 
properly normalized Lorentzian kernel tends to 8{ui— (Tr) and thus Lkk'icRjci/^ — » F kk i(an). 



In condensed matter calculations f8HlO| the Lanczos algorithm has been applied to the 
calculation of the Green function with a small value of 7], and its imaginary part has been 
interpreted as F kk i(uj) directly. This can be done if the spectrum is discrete (or discretized) 
and rj is sufficiently small. In our case we have a genuine continuum problem and we want to 
avoid any discretization, therefore we calculate L kk i{an) in the same way i.e. with finite ca- 
using the Lanczos algorithm, but then we antitransform L kk >(aR) in order to obtain F kk i{uj). 
Choosing 
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0|Vo> 



'(^olOtO^o) 

as starting vector for the Lanczos basis and setting z = Eq + <jr + ioi one finds 

L(a) = l^ o |6to|^ o )l m {(0 o |^^|0 o )} 
oi z - H 

showing that the LIT depends on the matrix element 

, , , 1 



(14) 



^00 = {<Po\ 



>■ 



z-H 

One can calculate xoo applying Cramer's rule to the solution of the linear system 

y^X z ~~ H) mn x n o = 5 m o 

n 

that arises from the expansion of the identity 

(z - H)(z - H)- 1 = I 



(15) 



(16) 



(17) 



(18) 



over the Lanczos basis {\<fii),i = 0, . . . ,n} PJTB|. In this way one is able to write a; o as a 
continued fraction containing the Lanczos coefficients and foj, 



^oo 



(z - Oq) 



(19) 



(z— oi)- 



(z-a 2 )-i)|.... 

and thus also the LIT becomes a function of the Lanczos coefficients 



L(a) = — Im < 
07 



(V>o|OtO|^ ) 



[z - Oq) 



> . 



{z-ai)- 



(z-a 2 )-trt 



(20) 



This illustrates that the Lanczos method allows to determine L(a) without solving the 
differential equation @, i.e. without inverting the hamiltonian matrix. 



III. EXCLUSIVE PROCESSES 

Different from an inclusive reaction where the final state of the studied particle system 
is unobserved, for an exclusive process one has to determine the cross section of a specific 
breakup channel for specific angles of the outgoing fragments. For instance, the unpolarized 
cross section for a two-body breakup induced by an external probe is given by 

-7^- = J2 c kk>(u,q,(fr f )f kk ,(uj,q,6 f ) , (21) 
dU f kk> 
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where f2/ = (Of, 0/) is the spherical angle of one of the outgoing fragments. In (^) the c kk * 
are known kinematical functions, while the dynamics of the particle system is contained in 
the structure functions 

fa{u>,q,O f ) = Y,T: km (E f )T ak , m (E f ) . (22) 

am 

The above transition matrix elements are defined as follows 

T akm (E f ) = (ipf, a \O k (q)\ipo, m ) , (23) 

where a stands for a set of quantum numbers of the final state function ipf with energy Ef 
in the specific breakup channel and m is the projection of the total angular momentum of 
the initial state wave function ipo. 

Different from the inclusive cross section, where an integral method leads to the direct 
determination of the dynamical function F kk /, one has to determine any single transition 



matrix element T akm (Ef) in order to calculate the exclusive cross section |TT . 

We illustrate the derivation of the integral transform method for the exclusive case 
restricting ourselves to a breakup of the A-body system in two fragments. In this case 
the final state of the system can be written using the Lippmann-Schwinger equation 

<fel = (€ W \ + (€ W \Vn 1 » , (24) 

Ef + ie — H 

where \<Pa W ) describes the plane wave relative motion of the two fragments with the same 
set of quantum numbers a as ipf, a , while V\% denotes the interaction of the constituents of 
one fragment with the constituents of the other fragment. In coordinate space the plane 
wave has the following form 

-ip-r 

€ W (r)=M^ 2 j^ , (25) 

where A is an antisymmetrization operator acting on the constituents of both fragments, ip\ 
and ip2 are the internal wave functions of the two fragments, while p is the relative momentum 
and r the relative coordinate of the two fragments. Note that internal coordinates and 
quantum numbers of ipi and ip2 are suppressed in (^) . 

Substituting ( |2~4] ) into ( |2"3"| ) one finds that the transition matrix element is described by 
the sum of two pieces, a Born term 

TSzrm = (€ w \Ok\^, m ) (26) 

and a final state interaction (FSI) term 

T£im = (€ W \Vi2- g; 6 k \^ m ). (27) 

Ef + ie — H 

While the Born term can be calculated without major problems the FSI term is more 
difficult. Using the completeness of the eigenstates \ip p (E p )) of H one gets for the FSI term 
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Tgi{Ej) = f df(^ w \V 12 \M E f)KM E r)\Tr^ ^-O^m) (28) 

J hif + ie — H 

df\^ w \Vu\Mm(MEf)\O k \^ m ) = - 1 = (29) 

hif -+- l€ — III f i 

-i7tF akm (Ef)+V f°° dE r F ; km{E J' ] (30) 



with 



F akm (Ef>) = T,(^\ V n\M E f')KMEr)\O k \ifa m ), (31) 
/' 

where the sum over /' stands for a possible degeneracy of the states \if)f) with energy Ep. 

Equation ( pi]) is only a formal solution for F ak m(E), since one should know all the 
eigenstates of H. A direct calculation of F akm {E) is even more difficult than the original 
problem of Eq. (f23|) , where one needs to know just the eigenstate of energy Ef. However, 
one obtains F akm indirectly applying the LIT approach: 

JEo-s (E — a)[E — a*) 

df{€ W \Vn^ -Ok\^o, m ) (33) 

Eo-S H — o H — o* 

= (€ W \Vu * * 6 h \^ m ) (34) 
H - a H - a* 

= (^aWl^tenW) • ( 35 ) 

The Lorentz state ipkmi.®) is obtained solving essentially the same differential equation as 
for the inclusive case discussed in the previous chapter: 

(H - a R + i<Ti)\$km) = O k \i^ 0>m ), (36) 

whereas V ; a( cr ) is the solution of new differential equation 

(H - a R + iai)$ a ) = Vn\€ W )- (37) 



If V\2 represents a finite range potential, the term on the right hand side of (|37| ) vanishes 
for large distances, and one has also for ip a an asymptotic bound state boundary condition 
similar to a bound state. If Vvz is not of finite range one can reformulate the LIT method 
accordingly [TJJjr5|]. 

Inverting the various L akm one gets the functions F akm and thus the transition matrix 
elements T£££. An additional calculation of the Born terms T^g™ leads to the searched 
structure functions. Also for exclusive cross sections one obtains reliable results with the 
LIT method as shown in a test calculation for the d(e, e'N) reaction [15 



Now we turn to the application of the Lanczos method to the LIT of the exclusive 
reaction. To this end we first rewrite the LIT in the following form 

L akm (a) = -L^IV^ 1 _ _ 1 —)d k \^ m ). (38) 

^cr i ii — <7r — tcrj H — <Jr + 1(7 1 
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Using the completeness relation over the Lanczos vectors {\4>i),i = 0, ...,n}, where the 
starting vector |0 O ) is the same as in (|14"D , and setting z = <Jr + iaj one finds 



i n 1 1 

Lakm(cr) = — {X!(0a W |^12|^t)[<0i| f-O k \lp 0m ) - (<^| T O k \lp ,m)] } 

2(7/ ~0 Z - H Z* - H 

= — ^ IZ^a Vl2 &} ^|0o - (0i — g: 0o) } 

^07 z - H z* - H 



i=0 



i\J {^,m\d ] k d k \^ m ) . 

iZ^Wa \V12mnxio -x i0 }}, (39) 



2(7 1 i=0 

with the matrix elements 

x 0i = <0o| |0i>, (40) 

5 i = (0o| -d&), ( 41 ) 

z* — H 

which can be written as continued fractions of the Lanczos coefficients, like in ([L9|). 

In the following we describe a simple algorithm for a recursive calculation of all the x$i 
starting from x o- Defining 

giy) = ^ (42) 

Z — a,j -2 — 



K+2 



one has 



z - a + g(l) 
The next matrix element Xqi can be written as 



^00 = 7TT • (43) 



XQ1 = 7 v\ — \ , x i r )\ ' ( 44 ) 

[z - ai)A - 61 + Ao</(2) 



where 



A = — — (45) 



is obtained from x o- Now a new parameter 



Ai = (2 - a f°- &1 (46) 



can be obtained from Xq\ leading to the evaluation of 



X ° 2 (z - a 2 )A x - Me, + A^(3) ■ (47) 



S 



Proceeding in this way it is easy to show that the n-th matrix element xo n can be written 



as 



with 



(z - a n )A n _i - n A n _2 + X n -ig{n + 1) 



^ _ {z — o n )A n _i — b n X n -2 
b n +i 



IV. RESULTS 

In order to check the accuracy and to quantify the advantages of the Lanczos method 
for the LIT approach we apply it to the calculation of total photoabsorption cross sections 
of 3 H and 4 He and compare to results which we obtain solving Eq. @ via an inversion of 
the hamiltonian matrix. 

For the total photoabsorption cross section we use the dipole approximation, where one 
has a single transition operator, namely 

O k = 6 v = 'tzJ-^- t (50) 
i 1 

here Z{ and r z (i) are the third components of the spatial coordinate and of the isospin of par- 
ticle i, respectively. We expand ground and Lorentz states on the hyperspherical harmonics 
basis making use of the recently developed EIHH approach [T^]. As NN interaction model 



we employ the semirealistic MTI-III potential ||17|| . We note in passing that calculations of 



the total photoabsorption cross section of 3 H and 4 He via the LIT already exist @J||7|. 

In Fig. 1 we show the LIT for the triton photodisintegration choosing the isospin channel 
T = 1/2 of the disintegrated three-body system. One sees that there is an excellent agree- 
ment between both calculations. The relative difference of both results is far below 0.1 %. 
For the second isospin channel (T = 3/2), not shown here, the comparison is of the same 
excellent quality. In Fig. 2 we show the LIT for the case of the 4 He photodisintegration, 
where one has only a single isospin channel of the disintegrated four-body system (T = 1) 
for the considered dipole transition. Again one sees a very good agreement comparing the 
results with and without application of the Lanczos method. One finds only tiny differences 
between both results. It is evident from Figs. 1 and 2 that the Lanczos method leads to 
highly reliable results. 

The main advantage of the Lanczos technique lies in the enormous reduction of the cpu 
time. In particular one should note that it is not necessary to perform all the Lanczos steps. 
In fact, checking the convergence of L(a) with respect to the number of Lanczos coefficients 
used in the continued fraction one finds the following: for <jr — *■ oo only the first Lanczos 
step is necessary, because the entire continued fraction is dominated just by gr. Also in 
the low <7r region one can expect that the continued fraction converges rapidly, since for 
or —>■ the dominating contribution to the term \j{z — H) is given by the lowest eigenvalue 
of H, which is well approximated using only the first few Lanczos coefficients. Thus one may 
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hope that in the intermediate <jr region a rather small number of Lanczos coefficients will 
be sufficient. In Fig. 3 we illustrate the convergence behaviour with respect to the number 
of Lanczos steps for the LIT of the 4 He photodisintegration for a test case with 1372 basis 
states. As one can see the transforms are somewhat different from the converged one only 
if one takes less than 100 Lanczos coefficients. The relative error is already smaller than 
0.5% for 150 coefficients and becomes rapidly smaller with a further increasing number of 
Lanczos steps. One can conclude that 200-^300 steps are completely sufficient in the studied 
case. Also due to the rapid convergence we obtain a reduction in cpu time of two orders of 
magnitude comparing with the conventional LIT calculation. 

Applying the Lanczos method one has to pay attention to the loss of orthogonality of the 
Lanczos vectors with an increasing number of steps, since it could reduce the precision of the 
Lanczos coefficients. We check the orthogonality calculating the scalar products of the 0-th 
and the i-th Lanczos vectors. The results for the LIT case, illustrated in Fig. 3, are shown 
in Table |. It is evident that after 300 Lanczos steps, where one has a sufficiently converged 
LIT, one has still a rather good orthogonality. In fact the overlap between the 0-th and 
the 300-th Lanczos vector is quite small (about 10~ 4 ). However, in other cases one could 
encounter problems with the orthogonality and thus one would need to reorthogonalize the 
Lanczos vectors. 

In the following we summarize our results briefly. The LIT approach for the calculation 
of both inclusive and exclusive cross sections including complete final state interactions is 
reformulated in order to use the Lanczos algorithm. The Lanczos method allows one to 
calculate the LIT without complete inversion of the hamiltonian matrix. It leads to a huge 
simplification of the calculation reducing the needed cpu time enormously (about two orders 
of magnitude) and, even more important, paving the way for calculations with large number 
of basis functions. We have tested the reliability of the method with great success for the LIT 
of the total photoabsorption cross sections of 3 H and 4 He. The high precision of the method 
together with the strong reduction of the cpu time opens the way to precise microscopic 
calculations of reaction cross sections also for more complex few-body systems. Indeed, very 
recently the new method has been applied to the first microscopic calculation of the total 
photoabsorption cross section of the A=6 nuclei with complete final state interactions of the 



six-nucleon system [IS 
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TABLES 



1 1.3 10" 16 

150 -2.4 10~ 6 

200 -1.2 10~ 4 

250 -9.5 10~ 5 

300 1.1 10~ 4 



TABLE I. Scalar product between the 0-th and z-th lanczos vector. 
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FIGURES 



FIG. 1. LIT of the 3 H photoabsorption cross section for T = 1/2 with conventional and Lanczos 
methods: absolute (a) and relative values (b) setting the conventional LIT result equal to 1 (07 = 20 
MeV). 

FIG. 2. LIT of the 4 He photoabsorption cross section with conventional and Lanczos methods: 
absolute (a) and relative values (b) setting the conventional LIT result equal to 1 (07 = 20 MeV). 

FIG. 3. Convergence of the LIT of Fig. 2 with respect to the number n of Lanczos coefficients: 
absolute (a) and relative values (b) setting the LIT result with 1372 Lanczos coefficients equal to 
1. 
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